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I summarize theory at Quark Matter 2002, stressing the continuing inabihty of a single 
^ . model to describe all notable features of the data from \/s/A : 55 200 GeV. 
> 

12 ; 1. INTRODUCTION 

o 

CN I The collisions of large nuclei may give us insight into the nature of QCD at high 
temperatures and/or densities. In this talk I review the status of theory for heavy ion 
O ' collisions from Quark Matter 2002. I concentrate on the the results for the largest nuclei, 
with atomic number A ^ 200. 

At the SPS at CERN, there are two notable results for AA collisions, for \/s/A : 5 
O : 17 GeV [|: 

J/\l' suppression: the number of J/\E' pairs is smaller in the most central collisions, 
versus the extrapolation from peripheral collisions, or from collisions with smaller A 0. 
The effect is most striking for the largest nuclei. 

Excess dileptons below the p: the rate of e"'"e~ pairs exceeds that in conventional 
\ hadronic models [0,0, although a broadened p-meson explains the excess. The effect 
is more prominent at lower, and not high energies, suggesting a density dependent effect. 
This also supports interest in going to even lower energies, such as at the proposed GSI 
collider. 

At BNL, RHIC has run at energies of ^/s/A = 55 GeV (briefly), at 130 GeV during 
Run I, and at 200 GeV during Run II. Results from Run I were first presented at Quark 
Matter 2001; those from Run II, at this Quark Matter, 2002. 

There is one notable change expected change between the SPS and RHIC. At the SPS, 
the particle multiplicity in AA collisions is a single peak about zero rapidity. By RHIC 
energies, a Central Plateau was expected to open up, in which physics is (approximately) 
boost invariant, independent of rapidity. Away from the incident nucleons of the fragmen- 
tation region, the Central Plateau is where a system at nonzero temperature, and almost 
zero quark density, might emerge. Perhaps even deconfined matter, as the Quark-Gluon 
Plasma. 

RHIC experiments find that the Central Plateau is rather narrow. In all, particles 
are spread out over k, ±5 units of rapidity. The multiplicites for identified particles are 
nearly constant over the central ~ ±1 unit of rapidity However, the pion's average 
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transverse momentum, pt, is only constant over ±.5 units of rapidity [||. But as with other 
variables, there is much to gain by looking at all rapidities. 

Even without (much) electromagnetic data, which proved so interesting at the SPS, so 
far there are four notable features of the RHIC data: 

"High"-pt suppression: the number of particles with transverse momentum pt between 
2 — 10 GeV is suppressed, relative to that in pp, times the number of binary collisions 
[3,P|, p!0tpT| , p!2[ . The overall suppression is by factors of 2 — 4. The suppression is now seen 
to be approximately constant for these pt 0. This is opposite what happens at the SPS, 
where high pt particles are not suppressed, but enhanced by factors of 2 — 3, through the 
Cronin effect. 

Elliptic flow: is a measure of momentum anistropy in non-central collisions. Hydrody- 
namics predicts elliptic flow is linear in pt for pions. This is seen up to ~ 1.5 GeV, as 
is the hydrodynamic behavior of protons. For : 2 — ^ 6 GeV, though, the (total) elliptic 
flow is flat [0,0, which is not expected from hydrodynamics (or anything else). 

HBT radii: pion interferometry gives a measure of the spatial size(s) of the system. 
Hydrodynamics predicts that a certain ratio of two sizes, Rout/Rside, is greater than one, 
and increases as pt does. Instead, experiment finds that Rout/Rside decreases with increas- 
ing Pi, and is about one by pt ^ 400 MeV [0,^. HBT radii indicate that hadronization 
occurs as a type of "blast" wave 111 



Jet absorption: At these energies, jets are seen in pp collisions, but an angular correla- 
tion finds that in AA collisions, the "backward" jet is strongly suppressed [jl6|. That is, 
in AA collisions there is stuff there which eats jets. 

It is thus an extremely exciting time. For now, experiment is triumphant: especially 
after the run at \fsl A = 200 GeV, there is a striking agreement between different experi- 
ments for many quantities of interest. The crucial advance has been precise measurements. 
If quantities are only measured to ±30%, then a lot of theories can hide under the error 
bars. With errors ±5%, it becomes possible to definitively rule (many) theories out. 

The challenge to theorists to synthesize all of these measurements into a consistent 
framework. Common prejudices held before the RHIC data — such as large increases in 
multiplicity, and large sizes from a strongly first order transition — are extinct. Certain 
models can explain some aspects of the RHIC data, but at present, no model can explain 
everything interesting. 

My review of theory at Quark Matter 2002 is organized thematically, emphasizing 
results of relevance to experiment. Consequently, I will offend by omission, neglecting all 
work on color superconductivity and other interesting topics. 



2. THE BEDROCK: THE LATTICE 

Our understanding of QCD at nonzero temperature rests upon numerical simulations 



on the Lattice . While dynamical quarks must be included in any realistic simulation 
of QCD, at present the Lattice is not near to the continuum limit with light, dynamical 
quarks. What it can do, with precision near the continuum limit, is tell us about the 
pure glue theory. While it is easy to bemoan the lack of nearly continuum data with 
quarks, it should not obscure this fundamental advance. That we can compute reliably 
the thermodynamics for non-abelian fields with three colors is no mean feat. 
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In the pure gauge theory, the deconfining phase transition is rigorously associated with 
the spontaneous breaking of a global Z{3) symmetry, above a transition temperature T^. 
Taking the string tension to be (400 MeV)^, ^ 270 MeV, within errors which are, at 
present, ~ ±5%. One expects a confined phase of hadrons below T^, and a deconfined 
Quark-Gluon Plasma above T^. The great surprise from the Lattice is that there is a third 
region, sandwiched between these two. This happens because the deconfining transition 
is so weakly first order that there is a nearly critical region about T^. For example, just 
below Tc, the string tension is about one-tenth its value at zero temperature. The nearly 
critical modes which become light about Tc are electric Z{3) glueballs, or Polyakov loops. 
Magnetic Z{3) glueballs, or 't Hooft loops, stay heavy. 

The Lattice also gives us insight into the transition as the number of colors is varied 
from three. It is truly second order for two colors, with a critical region whose width, in 
terms of the reduced temperature T/T^, is wider than for three colors. This leads me to 
the conjecture that for A*",, colors, while the transition is of second order when < oo, its 
width, in terms of T/Tc, shrinks like 1/Nc as Nc oo. At Nc = oo, then, the transition 
is of first order. 

For three colors, the nearly critical region is within ^ ±10% of T^. A consequence of a 
narrow critical region is that the potentials for Polyakov loops change very rapidly, within 
temperatures ~ ±2.5% of |jl8|. For the rapid expansion rates present in heavy ion 
collisions, this suggests that even if the system is thermal at Tc, its evolution below Tc is 
far out of equilibrium |]18 . 



A recent advance on the Lattice is to compute spectral functions using the Maximum 
Entropy Method [|1^ . The spectral density for the pion has a narrow peak in the confined 



phase, but is a broad peak in the deconfined phase, at least by 1.5Tc. In the p channel, at 
.6Tc the peak has moved up a little, and broadened more; the p is completely washed out by 
1.5Tc. Spectral functions for heavy quark bound states have also been computed. Notably, 
the Lattice finds that the J/"^ remains bound even at 1.5Tc, although less tightly bound 
states, such as the "if', evaporate below Tc. This is in contrast to potential models, which 
find, due to the decreasing string tension, that all charmonium bound states disassociate 
below Tc |0|. 

At present, the Lattice can only simulate dynamical quarks if the pions are too heavy. 
Estimates are Tc ~ 175 MeV, with no true phase transition in the thermodynamic limit, 
but only cross over behavior. It is possible that a first order transition reappears for 
physical pions, closer to the continuum limit [|17|. In principle, the transition can be more 
strongly first order than in the pure glue theory, since with the addition of three fiavors 
of massless quarks, the ideal gas pressure goes up by a factor of three. (It is useful to 
compare to the ideal gas pressure, since asymptotic freedom implies that the pressure is 
ideal at infinite temperature.) 

Present simulations with dynamical quarks demonstrate an approximate universality, 
termed "fiavor independence" . This is the observation that as a function of T/Tc, the ratio 
of the true pressure, to the corresponding value in an ideal gas of quarks and gluons, is 
nearly universal, independent of the number of fiavors. This is remarkable, given the large 
changes in both Tc and in the ideal gas term. This suggests that the thermodynamics of 
QCD, with three colors and three fiavors of dynamical quarks, is dominated by that of 
the pure glue theory. 
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Advances have been made with systems at nonzero (quark) chemical potential fi. The 
Lattice can generate a thermal distribution of quarks, but with standard Monte Carlo 
techniques, it cannot fill a Fermi sea. However, for a quark with energy E, if the temper- 
ature T is much greater than /i, then the Fermi-Dirac distribution function at 7^ is 
approximately that for /i = 0: 

1 1 

(1) 



That is, at high temperature it doesn't matter much whether or not you fill the Fermi 
sea. Thus the Lattice can compute at T fj,: by /i ~ 200 MeV, Tc has only decreased 



a small amount, to ~ 160 MeV This is well on the way to nuclear matter, which at 
zero temperature, occurs when /i exceeds fa 313 MeV. This is a new and strong constraint 
on effective models of thermal QCD, which often give too large a decrease. 

3. STATISTICAL MODELS: A BIG BOOST (VELOCITY) 

For the most central collisions at zero rapidity, an amazing summary of the single 
particle spectra is a thermal fit |l6|, |22| , ^ , p^ . Fits in which chemical freeze-out occurs at the 



same temperature as kinetic freeze-out are favored, with T ^ 165 MeV and /x ~ 14 MeV 



25|] . Resonance decays describe the excess of low momentum pions in central collisions. 
What happens for peripheral collisions, which must also have resonance decays, and yet 
do not have an excess of low momentum pions, is less clear, and usually ignored. 

The approximate equality of the temperatures for chemical and kinetic freeze-out is 
peculiar. Any scattering in a hadronic phase produces chemical freeze-out at a higher 
temperature than that for kinetic freeze-out p6| , p7i , so the data suggest that both tem- 
peratures are really one of hadronization, with little rescattering in a hadronic phase. 
This is one hint of possible non-equilibrium behavior at RHIC. 

The temperature for chemical freeze-out is consistent with data at lower y/s/A. From 
energies of y/s/A from a few GeV on up, chemical freeze-out occurs along a curve in which 
the energy per particle is constant, about ~ 1 GeV. In the plane of T and /x, even if the 
hadronization temperature agrees with Tc at /i = 0, it is distinctly lower than Tc{fi) for 
/i 7^ 0. For example, chemical freeze-out at AGS energies gives about fi = 200 MeV, and 
a hadronization temperature which is at most ~ 120 MeV; from the Lattice, though, at 
this n the corresponding is much higher, fti 160 MeV ET . 



To describe the behavior of particles with increasing mass, it is necessary to assume 
that all hadrons are emitted with respect to a local moving rest frame. At RHIC, the 
radial velocities of this local rest frame go up to ~ 2/3 c; averaged over radius, they are 
about ~ 1/2 c. This can be seen by eye: versus pt, single particle distributions for pions 
turn up, while those for protons (say) turn down. 

The radial dependence of the velocity of the local rest frame is not constrained by the 
data, and is fit to agree with the observed spectrum. The same is true of hydrodynamical 



models This is why they are fits. For example, consider how the single particle 

spectra change with rapidity, or centrality. While the temperature might be the same, 
the local fiow velocity now depends not just upon the radius, but also upon the rapidity, 
centrality, etc. It is untenable to consider only zero rapidity, and ignore the rest. 
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A statistical model implies not only what the chemical composition is, but, as well, the 
Pt-dependence of the single particle spectrum. Of course a thermal distribution should 
only hold up to some upper scale, perhaps 1 — 2 GeV. It would be interesting to compute 
the ratios of moments of transverse momenta: 



th 



(2) 



Here exp and th denote, respectively, moments computed from experiment, versus a 
thermal distribution (with some assumed velocity profile). By definition, if the overall 
number of particles is thermal, ro = 0. For n > 1, r„ is a dimensionless series of pure 
numbers; the fit is good until r„ is no longer small. This must happen at some large 
n, since eventually fluctuations from hard momentum processes dominate. It would be 
interesting to determine these ratios from experiment, for all collisions in which a thermal 
fit works. 



4. "IDEAL" HYDRODYNAMICS AND ELLIPTICAL FLOW 

A dynamical realization of a thermal fit is a hydrodynamical model. A specific equation 
of state is assumed, with some initial temperature and velocity profile at some initial time. 
The system then evolves according to an equation of state, with the parameters chosen 
to agree with the observed single particle spectra. First order transitions are usually 
assumed, but this doesn't matter much, given the rapid expansion characteristic of heavy 
ion collisions. 

A measure of hydrodynamic behavior is given by elliptic flow. For a peripheral collision, 
in which the two nuclei only partially overlap, an "almond" is formed in the plane per- 
pendicular to the reaction plane. As the system hadronizes, this spatial anistropy turns 
into a momentum anistropy, with the average momentum larger along the narrow part of 
the almond then along the long part. This elliptical anistropy |I^p!3| has been measured 



as a function of centrality and pt'-, overall, the values at RHIC are about twice as large as 
at the SPS. By geometry, elliptic flow vanishes for zero centrality, as nuclei which com- 
pletely overlap cannot have any anistropy. Elliptic flow, as measured from two-particle 
correlations, fails to do this, so that more sophisticated measures, involving correlations 
between four or more particles, are imperative . 



Hydrodynamic models predict that for pions, the elliptic flow depends linearly on the 
transverse momentum. The local flow velocity also predicts the behavior of elliptic flow 
for heavier particles, such as protons. Both predictions are borne out by the experimental 
data, for momenta up to pt ~ 1.5 GeV. 

Versus centrality, as measured by the number of participants, hydrodynamics predicts 
that the elliptic anistropy is linear near zero centrality, which is observed. When the 
number of participants is half the maximum value, though, hydrodynamics significantly 
overpredicts the elliptic flow. 

Calculations with three dimensional hydrodynamics, giving single particle distributions 
versus rapidity, have been performed [^. Even with agreement at zero rapidity, the 
results from these hydrodynamic calculations are much broader in rapidity than the ex- 
perimental data [ll3| , p!^ . 
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There are two conundrums about these hydro dynamic fits: first, the initial times re- 
quired in hydrodynamic calculations are extremely small, ~ .6 fm/c [^. Secondly, they 
assume ideal hydrodynamics. Even with times expected from saturation (see below), this 
is a very short time. Further, in QCD viscosity coefficients can be computed from the 
Boltzman-Altarelli-Parisi equations. At least in weak coupling, the shear viscosity is large 
(on a natural scale), and cannot be neglected |T^J^. 

These technicalities should not obscure the fact that elliptic flow demonstrates that 
heavy ion collisions exhibit significant collective behavior. Further, while fits to single 
particle spectra with ideal hydrodynamics and short times may well work, this does not 
exclude the possibility of similar fits with non-ideal hydrodynamics; it is not clear how 
the times might then change. 



As an example, [0 assumes that particles obey a thermal distribution in the transverse 



direction, but free stream in the longitudinal (or beam) direction. Adjusting the initial 
conditions to give the right single particle distributions, the resulting elliptical anistropy 
is too small by a factor of two. As in many other examples, the data sharply limits theory. 

Experimentally, it is unremarkable that hydrodynamics fails above pt ^ 1.5 GeV. Hy- 
drodynamics should break down at short distances; that it works down to ~ .13 fm is 
actually pretty good. Rather, the surprise is that the elliptic anistropy is approximately 
constant for pt : 2 —>■ 6 GeV. In QCD, one expects cross sections to peak at some momen- 
tum scale on the order of a few GeV, and then to fall off with the powers characteristic 
of QCD. It is very difficult to imagine how anything fiat in pt could ever emerge. 

One possibility is that one is not measuring collective phenomenon at all, but simply 
some property of two jets, which is misinterpreted as collective flow. used a mini-jet 
model to obtain an elliptic anistropy which agrees with experiment; the single-particle 
spectra, though, is wrong. Methods to determine elliptic anistropy using correlations 
between four or more particles should be insensitive to contamination by jets . 



5. HBT RADII: A "BLAST" WAVE? 

For identical particles, a length scale can be determined by pion interferometry through 



the Hanbury-Brown-Twiss (HBT) effect |T5|]. This length scale is related to the surface 



at which the pions last interacted. Since there is axial symmetry to a heavy ion collision, 
there are three distances, corresponding to along the beam direction, Riong, along the hne 
of sight. Rout, and perpendicular to that, Rside- 

One of the big surprises from RHIC is that the HBT radii did not grow much between 
y/s/A = 17 to 200 GeV. The change in RiongRsideRout is, more or less, the same as the 
increase in multiplicity, ~ 50%. 

This can be taken as direct experimental evidence for the absence of a strongly first 
order phase transition in QCD, completely independent from the Lattice. If the transition 
were strongly first order, as it went through Tc the system would supercool and grow in 
size. Estimates of the sizes of the system before QM'Ol ranged up to tens of fermi, which 
are not seen. Unfortunately, putting a bound on the latent heat of the transition is 
manifestly a model dependent exercise. Still, it would be an amusing exercise. 

The details of the HBT radii, however, have proven to be much more interesting than 
expected. Before the RHIC data, it was thought that the hadronic firetube from an 
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AA collision might be like a "burning log" . But instead of smouldering, the RHIC data 
suggests that the log blows up. 

In particular, the results from RHIC appear to contradict any hydrodynamic description 
|]I5|,^. Versus experiment, hydrodynamics gives values of Riong and Rout which are 
too large, and a Rside which is too small. Of especial interst is the ratio of Rout/Rside- 
hydrodynamics predicts this ratio should be ~ 1.5 2, and which increases with pt. At 
RHIC, the ratio decreases as pt goes up, and is about one, .85 > Rout/Rside > 1-15 [|15]- 

The HBT data can be parametrized as a type of "blast" wave, with a velocity ~ 3/4 c 
jl5|. This may indicate a type of "explosive" behavior |T8| , |33[| , a term first used by the 



experiment alist s . 

HBT radii can also be used to give a measure of the entropy carried by pions |I5| . 
Per degree of freedom, the entropy in peripheral collisions is greater than that in a Bose- 
Einstein gas, but is less than that for the most central collisions. This may indicate 
coherence in the initial state. 

Partonic models in which the cross section is artificially multiplied by a factor of ~ 14 



34| give both Rout/Rside ~ 1, at least by pt ~ 1.5 GeV, and a fiat elliptical anistropy. 



One may view this as evidence that at hadronization, apparently near T^, the partons 
appear to behave much more strongly than expected from perturbative QCD. 



6. "HIGH"-j9i SUPPRESSION 

From the first RHIC data, it was clear that the spectra for "high"-pf particles, meaning 
above, say, 2 GeV, is qualitatively different in central AA collisions, versus pp collisions 
at the same energy. Dividing by the number of participants, the number of particles 
at high-pf is significantly less in central AA collisions than in pp, by overall factors of 
2 — 4 [P,p|, p^pT|p]^ . This is quantified through the ratio Raa, which is the ratio of the 
number of particles in central AA collisions, divided by that in pp, as a function of pt- 
The suppression begins above pt ^ 2 GeV; above 4 GeV, Raa ~ 1/3 — *■ 1/4 for charged 
hadrons, and Raa ~ 1/5 — >■ 1/6 for pions [§. A surprise of the Run II data is that for 

: 2 — i> 9 GeV, Raa is approximately constant, up to at least 9 GeV 0. 

This suppression of high-p^ particles is opposite to what is observed at the SPS. There, 
due to what is known as the Cronin effect, the ratio Raa is greater than one, going up 
to ~ 2.5 by P( ~ 3 GeV. This change in the spectrum must be considered as one of the 
most dramatic features of the RHIC data. 

The usual explanation of high-p^ suppression is due to energy loss |lT0| , pJ] , p!2| . Bjorken 



originally noted that a fast quark (or gluon) loses energy as it traverses a thermal bath, in 
just that same way that any charged particle does in matter. Single particle distributions 
can be explained using parton models jl^ . 

Detailed features of the data appear difficult to explain as energy loss, though. Naively, 
one might expect soft processes to scale as the number of participants, ~ A, and hard 
processes to scale as the number of collisions, ~ A'^^^. As a hard process, then, one would 
expect energy loss to scale as A'^^^. Instead, in peripheral collisions, Raa scales with the 
number of participants [^. It is possible that the difference of two terms of order A'^^^ 
is of order A, but that seems unnatural. 

Further, the observed constancy of Raa for : 2 — 9 GeV is surprising; perturbative 
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models of QCD do not give constant behavior. The apparent constancy also reflects 
changes in particle composition, while pions dominate below pt ^ 2 GeV, unlike pp 
collisions, there are as many protons as pions above pt ^ 2 GeV. Taking protons as color 
baryon junctions, in perturbative QCD the pion part of Raa peaks at a smaller pt than 
the proton part p6 |. 

In weak coupling, energy loss in QCD is proportional to 

m 



of the jet 



E, where E is the energy 
If one assumes, simply for the sake of understanding the data, that the 
energy loss is proportional instead to E, with hard particles losing ^ 7% of their energy 
per scattering, then a next-to-leading order calculation in QCD gives a good fit to the 
data 137 . 



7. SATURATION 

Another surprise from the first RHIC data was that the multiplicity did not grow as 
rapidly as predicted, at least on the basis of various cascade models. A natural explanation 
for this is given by models of saturation [[Tll ,^. 39 , 40|. 41 . 42 1 . 



The application of saturation to AA collisions is, at the most basic level, purely a 
kinematical effect. Consider a nucleus-nucleus collision, in the rest frame of one of the 
nuclei. For atomic number A ^ 200, in its rest frame the incident nucleus has a diameter 
no greater than 2A^^^ 15 fm. By Lorentz contraction, this distance gets shrunk down 
by a factor which is about 1/(^5/^4). Eventually, the color charge of the incident nucleus 
looks not like a nucleus, but just like a very thin pancake, with a big color charge ~ A^^^. 
Assuming that distances on the order ofl/3— >l/4fm are small on hadronic scales, the 
incident nucleus looks like a thin pancake when y/s/A : 45 ^ 60 GeV. It is amusing that a 
simple estimate gives an energy right near where a Central Plateau, in which the particle 
density is constant with rapidity, first appears. 

In detail, saturation is a dynamic criterion. It states that at sufficiently small Bjorken- 
X, quark and gluon distribution functions are dominated by gluons, which peak at a 
characteristic momentum scale, termed the "saturation" momentum, Psat- (This gluon 
dominance is reminiscent of flavor independence for thermodynamics.) For any pertur- 
bative approach to work, psat cannot be less than at least 1 GeV. The above kinematic 
argument suggests that Pg^^i ~ A^/^: thus one can probe smaller x values with large nuclei 
at RHIC, say, than in ep collisions at HERA. 

What is most important about saturation is, again, almost a kinematical effect: it resets 
the "clock" for heavy ion collisions. In the Bjorken picture which dominated before the 
RHIC data, one assumed that hadronization occured at time scales ~ 1 fm/c; after all, 
what other time scale is there? Thus in the Bjorken picture, there seemed as if there was 
little time for even the largest nuclei, only 7 fm in radius, to thermalize. (Unless, again, 
there were a strongly first order transition, which is why it was so popular before RHIC.) 

With saturation, however, the natural scale of the clock is given by l/psat', for psat ~ 
1 GeV, this is already ~ .2 fm/c; see, e.g., [^. That is, saturation makes the hadronic 



"clock" runs at least five times faster! The possibility of interesting things happening is 
far more likely. 



Saturation is realized in the Color Glass Model |11,38,39|. The gluon fields from the 



incident nucleus are described as classical color sources, reacting much quicker than the 
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fields in the target nucleus. Taking a gluon field to scale with the QCD coupling constant 
(yf as A° ~ 1/(7, one concludes that the action, and indeed all quantities — such as particle 
multiplicity, average energy, etc. — scale like l/g"^. In an asymptotically free regime, 
then, all quantities grow like ~ l/as(Psat) ~ log(psat)- This small, logarithmic growth in 
the multiplicity agrees qualitatively with the RHIC data (although one really needs the 
increase to LHC energies to make this quantitative). 

This picture is only approximate. Even if a gauge field is ~ l/f/, the action need not 
scale like 1/ g"^. In AA collisions, at initial stages there is a screening mass generated along 
the beam direction, but not transverse, with a mass squared ~ as at leading order |lTO| , pJ| . 
Such a dynamically generated mass scale changes integral powers of l/a^ ~ \og{j)sat) into 
fractional powers; see, e.g.^ [p!0|JTl[] . 

Modulo these theoretical quibbles, it seems plausible that saturation describes the initial 
state of AA collisions at high energies. Fits to the particle multiplicity, including the 



dependence upon centrality and rapidity, agree approximately with the data [38|. It is 



not evident how to turn gluons into hadrons, as sometimes the mysteries of "parton- 
hadron duality" are invoked. It is surprising that such models work over a wide region of 
centrality and rapidity, since saturation (valid at small Bjorken-x) should not work well 
in the fragmentation region (which is large Bjorken-x). The particle density in such fits 
has peculiar peak at zero rapidity . Experimentally, the fragmentation region appears 
to be described by a universal limiting form ||3^ (although it is crucial to have data with 
identified particles, which is presently lacking). 

Saturation does not describe other basic features of the data, though. The most serious 
problem is that the average pt in saturation is {pt) ~ '^Vsat] even with psat as low as 1 GeV, 
this is an average pt ^ GeV. In contrast, at RHIC {pt) ~ 550 MeV. The average energy 
from saturation will decrease due to inelastic processes and the generation of entropy. 
Assuming that this fixes the overall constant, one is still at a loss to explain why the 
average pt changes by at most 1% between \/s/A = 130 GeV and 200 GeV, while the 
multiplicity changes by at least 15% 0. In saturation models, the average pt grows with 
multiplicity. 

A related problem is the chemical composition. Parton-hadron duality is really gluon- 
pion duality; but if the average gluon momentum is large, why don't the hard gluons 
become kaons? Instead, at RHIC kaons are much less numerous than pions, only about 
15% as much. 



The elliptic anisotropy can be calculated in a Color Glass [40|. As a function of mo- 



mentum, it is dominated by very soft momentum, peaking at a value of pt ~ psat/4:. This 
is not seen experimentally. 

Of course if saturation describes the initial state, and not the final state, then there is no 
problem with the above features of the data. Estimates of thermalization in perturbative 
QCD give numbers ~ 3 fm/c 0)1101 • Further, thermalization occurs not due to elastic 
processes, but is dominated by inelastic processes P3| . 

It is claimed that the logarithmic growth of the total cross section follows from satu- 
ration . This is a matter of black and white: if there are massless particles, there is 



no bound on the total cross section. Thus saturation, which has no mass gap, cannot 
guarantee a finite total cross section. 
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